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We estimate the radiative energy-loss of heavy quarks, pro- 
duced from the initial fusion of partons, while propagating 
in a quark-gluon plasma which may be formed in the wake 
of relativistic heavy ion collisions. We find that the radia- 
tive energy-loss for heavy quarks is larger than the collisional 
energy-loss for all energies. We point out the consequences 
on possible signals of the quark-gluon plasma. 



One of the most interesting predictions of QCD is 
the transition from the confincd/chirally broken phase 
to the deconfmed/chirally symmetric state of quasi-free 
quarks and gluons, the so-called quark-gluon plasma 
(QGP). Relativistic heavy ion collisions are being studied 
with the intention of investigating the properties of the 
QGP [Q. While experiments at AGS and SPS continue, 
new experiments have been planned at RHIC and LHC 
with centre of mass energies 200 AGeV and 5.5 ATcV, re- 
spectively. During the past decade many different signa- 
tures of the transition to the QGP have been proposed. A 
promising example is the emission of penetrating probes 
such as dileptons and single photons, which can reveal 
the early parton dynamics and the history of evolution 
of the plasma. Similarly, the production and propagation 
of open charm and high energy jets in a dense medium, 
can provide information Q about parton scattering and 
thermalization of the partonic system. Jets are expected 
to show up at collider energies at RHIC and LHC. 

Heavy quark pairs are mostly produced from the initial 
fusion of partons (mostly from gg — > QQ, but also from 
qq — > QQ, where q denotes one of the lighter quarks and 
Q is a heavy quark) of the colliding nucleons and also 
from the QGP, if the initial temperature is high enough, 
which is likely to be achieved at RHIC and LHC en- 
ergies. The charm quarks will be produced on a time 
scale of l/2m c ~ 0.07 fm/c, which would be as low as ~ 
0.02 fm/c for bottom quarks. There is no production of 
heavy quarks at late times in the QGP and none in the 
hadronic matter. Thus, the total number of heavy quarks 
gets frozen very early in the history of the collision which 
makes them a good candidate for a probe of the QGP. 
Immediately upon their production, these heavy quarks 
will propagate through the deconfined matter and start 



losing energy. The energy-loss suffered by these quarks 
will determine the shape of the dilepton spectra produced 
from correlated charm (or bottom) decay which provides 
a large background to dilepton production from annihi- 
lation of quarks in the plasma. We shall come back to 
this aspect towards the end of this Letter. 

There are two contributions to the energy-loss of a 
heavy quark in the QGP: one caused by elastic collisions 
with the light partons of the QGP and the other by radia- 
tion of the decelerated color charge, i.e., bremsstrahlung 
of gluons. There is an extensive body of literature PHq] 
on the collisional energy-loss of energetic quarks con- 
sidering clastic collisions with the quarks and gluons 
(Qg — > Qg and Qq — > Qq) of the dense medium. A com- 
plete leading order result for the collisional energy-loss 
of heavy quarks has been found using the hard thermal 
loop resummation technique || . 

It is well known that the contribution of the radia- 
tive processes (Qq — ► Qqg and Qg — > Qgg) is of the 
same order in the coupling constant as the collisional 
energy-loss §. The estimate of the radiative energy- 
loss in the past has been discussed by a number of au- 
thors 10 -13 within perturbative QCD taking into ac- 
count the Landau-Pomeranchuk suppression due to mul- 
tiple collisions. These studies, however, were limited to 
the case of massless energetic quarks and gluons. As far 
as we know, there is no estimate of the radiative energy- 
loss for heavy quarks in the literature. It is not easy to 
extend the sophisticated treatment of multiple scatter- 
ing formulated by the authors of Ref. p| to the case of 
heavy quarks. 

Until such a detailed investigation is performed, an ex- 
tension of the work of Ref. [M for the radiative energy- 
loss of massless quarks to the case of heavy quarks can 
provide valuable insight. This approach is similar to the 
one by Gyulassy, Wang and Pliimer p2] and leads to al- 
most identical results in the case of light partons. Baier et 
al. |l3| , on the other hand, found a different dependence 
of the radiative energy-loss on the energy of the parton 
by including the rescattering of the emitted gluon in the 
QGP. However, inserting typical values for the energy of 
the parton, the temperature, and the coupling constant 
yields quantitatively similar results. Therefore we will 
restrict ourselves to the simple approach of Ref. EG] in 
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the present work and discuss the consequences of our es- 
timate for signatures of the QGP. 
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FIG. 1. Feynmann diagrams for gluon bremsstrahlung 
from quarks. 

We start from an expression for the gluon emis- 
sion probability which has been derived by Gunion and 
Bertsch E3] in the case of light partons assuming a fac- 
torization of the matrix elements of Fig. 1 into elastic 
scattering and gluon emission. It can be shown that their 
result also holds if one of the light quarks in Fig. 1 is re- 
placed by a heavy quark assuming that the engery of the 
emitted gluon is not too large (qo <C y/s). The final re- 
sult for the multiplicity distribution of the radiated gluon 
can then be written as 
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where q = (qo, q±, q^) and I = (Z , ^3) are the four 
momenta of the emitted and the exchanged gluons, re- 
spectively, and r\ = (l/2)ln[(g + <?3)/(<7o - 93)] is the 
rapidity. C'a — 3 is the Casimir invariant of the adjoint 
representation. The factorization in ([!]) was obtained in 
the limit xlj_ << qj_, where x is the fractional momen- 
tum carried by the radiated gluon relative to the max- 
imum available, such that the radiation is confined to 
a uniform (central) rapidity region, the most important 
zone in relativistic heavy ion collisions. Eq.(|l|) holds also 
for gluon bremsstrahlung emitted by a gluon. 

Now we can estimate the radiative energy-loss per unit 
length for heavy quarks by multiplying the interaction 
rate T and the average energy-loss per collision v. which 



is given by the average of the probability of radiating 
a gluon times the energy of the gluon. One can further 
correct for the Landau-Pomeranchuk effect by including a 
formation time restriction [ p^|JTl| ] through a step function 
6(t — Tf). This puts a restriction on the phase space 
of the emitted gluons in which the formation time, rt 
must be smaller than the interaction time, r = 1/r. The 
formation time is estimated by requiring the separation 
between the emitted gluon and the parton from which it 
is emitted to be r± = v±t > l/q± (q± = \q±\) according 
to the uncertainty principle. Using v± = q±/qo and go = 
q±coshr], we find Tf = coshr]/q±. 

The average radiative energy-loss per collision is cal- 
culated as 



v = ( n g q ) = / dr/d 2 q ± 
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Performing the integration in (Q) in the limit (q±r) 2 >> 
1 and q± >> Zj_, we get 
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For the infrared cut-off g™ m we choose the Debye screen- 
ing mass of a pure gluon gas, 



qT n = VD = V^oTsT , 



(4) 



where T is the temperature of the system. For heavy 
quarks of mass M, the square of the maximum transverse 
momentum of the emitted gluon is given by 
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where s is the Mandelstam variable. To evaluate (|^) we 
need to compute (s) and (1/s) leading to 
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where E and p are the energy and momentum of a in- 
coming heavy quark and p' is the average momentum of 
the light quark or gluon of the QGP. The average value 
of p' can be taken as ~ 3T. Now, (||) becomes, 
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The average momentum transfer of the scattering process 
is defined as 
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where the differential cross section for elastic scattering 
is 
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and the square of the matrix element, \M.\ 2 can be ob- 
tained from Ref. p5[ . In the limit \t\ << s, the differen- 
tial cross section in (0) can be approximated by 
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We have checked that the modification of the energy- 
loss using the full expression for \Ai\ 2 is negligible (see 
below). Combining (JsJ) to ( [To| ) we get 
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The radiative energy-loss for heavy quarks is then ob- 
tained by combining (||) and (|ll|) and multiplying by 

r = l/r, 
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Since the mass of the quark in this expressions enters 
only via the maximum transverse momentum (Q) the ra- 
diative energy-loss of a heavy quark differs from the one 
of a massless quark only for small energies of the order 
of M. 

Let us also recall the expression for the collisional 
energy-loss of heavy quarks considered in Ref. j| using 
the hard thermal loop resummation technique. In the 
domain E « M 2 /T, it reads 



dE\ 
dx A 



inazT 



1 



6 



x In 





In 









1 



2 6+Tl fB(v) 



2v 2 
ET 



m g M 



(13) 



whereas for E » M 2 /T, it is 
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where v is the velocity of the heavy quarks, B(v) is a 
smooth function of v, which can be taken approximately 
as 0.7, rif is the number of light quark flavours taken 
as 2.5, and m g — •v/(I + nj /6) /3g s T the thermal gluon 
mass. 

It should be noted that the collisional and the radia- 
tive energy-loss are of the same order in the coupling 



constant although the latter is caused by higher or- 
der diagrams within naive perturbation theory. The rea- 
son for this behaviour is the fact that the interaction 
rate entering into the radiative energy-loss suffers from a 
quadratically infrared singularity using a bare propaga- 
tor for the exchanged gluon, whereas this divergence is 
reduced to a logarithmic one for the collisional energy- 
loss. This reduction is caused by the presence of the 
energy transfer of the exchanged gluon in the definition 
of the collisional energy-loss [[|. In the case of the ra- 
diative energy-loss, on the other hand, this factor is ab- 
sent, because the energy-loss is caused by the emitted 
and not by the exchanged gluon. Using a hard thermal 
loop resummed propagator the quadratic singularity in 
the interaction rate is reduced to a logarithmic one and 
the final result is of higher order (r ~ a s ) than naively 
expected. Multiplying this rate by the gluon emission 
probability (|l|) yields the result (|l^) of order a 2 . Using a 
resummed propagator for the collisional energy-loss leads 
to the finite expressions (|l3|) and ( |l4|) of the same order. 
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FIG. 2. The threshold momenta of the heavy quark as a 
function of T. 

The expression for the radiative energy-loss in (|l^) ex- 
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hibits a threshold behaviour: for (<7™ ax ) 2 < n 2 D there is no 
radiative energy-loss which is shown by the hatched area 
in Fig. 2. (Obviously, this behaviour will be different for 
mass-less quarks.) We see that the value of the threshold 
momenta of the heavy quarks, below which there is no 
radiation, increases with increasing temperature. 

In Fig. 3, we compare our results with that of the colli- 
sional energy-loss for heavy quarks obtained in Ref. || as 
a function of energy at a temperature T — 500 MeV and 
a s =0.3 for charm and bottom quarks. Before discussing 
our results, we give the justification for the approxima- 
tion made in ( |l0| ) for the differential cross section used in 
computing the radiative energy-loss, which enabled us to 
obtain the closed form given above. The solid lines rep- 
resent the radiative energy-loss of heavy quarks with full 
\Ai\ 2 whereas the dashed lines correspond to that with 
the approximate expression. We see that retaining only 
the ~ 1/t 2 term in da/dt is sufficient for our purpose. 
The dash-dotted lines in Fig. 3 represent the collisional 
energy-loss obtained in Ref. ||. We find that the radia- 
tive energy-loss dominates over the collisional one at all 
energies. For E > 20 GeV the difference amounts to an 
order of magnitude. 
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The QGP expected to be produced at RHIC and LHC 
is likely to be far from chemical equilibrium, initially. 
Chemical reactions among the partons will then push it 
towards a chemical equilibrium jl6]-[l8| . The evolution of 
the temperature and the quark and gluon fugacities r] at 
RHIC and LHC eneriges has recently been obtained p^] 
for such a scenario with initial conditions from a Self 
Screened Parton Cascade Model |l9[] . As a first esti- 
mate, the expressions for the collisional energy-loss given 
earlier (Eqs. ([l3|), fll4|)) can be modified by replacing the 
terms (1 + n//6) by (\ g + A g n//6) and n//(6 + nj) by 
\ q rif/(6\g + XqTif) to account for the departure from 
chemical equilibrium. Alternatively one may use the re- 
sults of Ref. pQ ] for a non-equilibrium plasma. The radia- 
tive energy-loss is modified by using the non-equilibrium 
Debye mass Il7j] 



4Tr\ g a s T 2 . 



(15) 



in (yp. 

This has interesting consequences. It has been shown 
recently |2(i(| that considering only the collisional energy- 
loss, in this manner, amounts to having only a small drag 
on the motion of heavy quarks in such a plasma, at least 
at RHIC energies, where the charm quarks were found to 
loose only ~ 10% of their energy during their propaga- 
tion, through the plasma and upto 40% of their initial en- 
ergy at LHC, in a collision involving two gold nuclei 

The drag acting on the heavy quark is conveniently 
defined by writing 



dE 
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where A denotes a drag-coefficient in the spirit of the 
treatment used earlier in literature (J^^] and p is the 
momentum of the heavy quark. Adding the collisional 
and radiative energy-loss experienced by a heavy quark in 
such an equilibrating and cooling plasma, we have found 
that 
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FIG. 3. The energy-loss of heavy quarks as a function of 
their energy. 



A ~ C/t, 



(17) 



where C is a slowly varying function of p with C ~ 0.4 for 
charm quarks at RHIC energies and ~ 0.7 at LHC ener- 
gies, for E < 5 - 6 GeV. This leads to a rather large drag 
of - 1.6/fm at RHIC and - 2.7/fm at LHC on charm 
quarks at r = r,, where Tj = 0.25 fm /c is determined by 
the onset of the kinetic equilibrium ]18| . 

This has a very important implication. Consider a 
charm quark having an energy Ei of the order of a few 
GeV at time Tj in such an expanding and chemically equi- 
librating plasma. Due to this large drag, the charm quark 



x The fugacities are defined |17| as /, = \ifi, where / is 
the distribution of the partons and / is the corresponding 
equilibrium distribution. 
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produced initially will come to rest very quickly and dif- 
fuse. We have verified this by performing numerical cal- 
culations of the final momentum of charm quarks which 
propagate under such a drag. We find that, irrespectively 
of the initial energy (Ei < 5-6 GeV) the final energy of 
the charm quark is about 1.5 - 1.6 GeV, in the cases 
considered here. Recall again that the charm quarks do 
not come to a stop if only the collisional energy-loss is 
included @. 

Thus we conclude that the radiative energy-loss of 
heavy quarks produced initially in relavisitic nuclear col- 
lisions plays a dominant role in pulling them to a stop 
in the QGP both at RHIC and LHC energies. Their fi- 
nal momentum distribution will then be determined by 
the temperature at which the hadronization takes place. 
This could be the temperature of the mixed phase, if such 
a phenomenon takes place. 

We may add that Svetitsky and co-workers [|| have 
actively investigated such a scenario. In their work the 
large drag coefficient arises due to a large value of a s ~ 
0.6 and a fully equilibrated plasma, even though only the 
collisional energy-loss is included. 

Let us return to the discussion of the momentum dis- 
tribution of charm quarks. (Similar considerations hold 
for bottom quarks.) It is expected that the momentum 
distribution of charm quarks will be reflected in the mo- 
mentum distribution of charmed mesons, whose corre- 
lated decay will provide a back-ground to dileptons from 
quark annihilation. We see immediately that a look at 
the p± distribution of these leptons may help us to isolate 
the two contributions, as they should be very different for 
the two sources. Shuryak and Lin et al. j?3| have ar- 
gued that the correlated charm decay back-ground for 
dileptons may be suppressed if the energy-loss of charm 
quarks is taken as 1-2 GeV/fm. Our study lends a strong 
support to their conclusion which were obtained by at- 
tributing an arbitrarily assumed value for the energy-loss. 

In conclusion, we have estimated the radiative energy- 
loss of heavy quarks propagating in a quark gluon plasma. 
This, along with the (fairly small) collisional energy- 
loss acts as a strong drag force on heavy quarks, which 
pulls them to a stop even in a chemically equilibrating 
and cooling plasma. This ensures that the momenta of 
the resulting charm mesons will be determined by the 
hadronization temperature. The correlated decay of such 
charm mesons will then no longer pose a back-ground for 
dileptons having their origin in the quark-antiquark anni- 
hilation at least at large invariant mass. This separation 
could even be made easier by measuring the p± distrib- 
tuion of the lepton pairs. In a future publication we shall 
report the result of the transverse hydrodynamic flow of 
the plasma on these conclusions [[[ 
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